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Abstract 

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy 
survey over a blank field as large as 34' x 27'. The program consists of very deep multi-band optical imaging, 
near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific 
goals of the project are to construct large samples of Lyman-brcak galaxies at z ~ 4 — 5 and Lyman alpha 
emitters at z ~ 5.7 and 6.6, and to make detailed studies these very high-rcdshift galaxy populations. 

In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked 
images and object catalogs in seven bandpasses. The optical imaging was made through five broad-band 
filters, B, V, R, i' , z' , and two narrow-band filters, NB816 (Ac = 8150A) and NB921 (Ac = 9196A) with 
almost 10 hours long integrations for each band. The limiting magnitudes measured at 3cr on a 2" aperture 
are B ^ 28.45, V = 27.74, R = 27.80, i' = 27.43, z' = 26.62, NB816 = 26.63, and NB921 = 26.54 in the AB 
system. The object catalog constructed for each of the seven bands contains more than 10^ objects. The 
galaxy number counts corrected for detection incompleteness and star count contribution are found to be 
consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have 
been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/. 
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1. Introduction 

Deeper and wider-field observations are required to un- 
derstand the universe in longer baseline of time and space. 
A large galaxy survey is a simple and essential method in 
today's astronomy for the purpose. However, an unprece- 
dentedly large galaxy survey is made possible only with 
the combination of some distinctive resources, such as a 
large telescope, a unique and efficient instrument based 
on advanced technologies, and a large amount of telescope 
time created by the copperation of many people who ap- 
preciate the value of the survey. In particular, such an 
extremely deep and wide galaxy surveys to collect a sta- 
tistically robust number of high-z galaxies beyond z = 3 
are definitely needed. There are a lot of high-z galaxy 
surveys so far, for Lyman break galaxies (e.g., Steidel et 
al. 1996; Steidel et al. 1999, Iwata et al. 2003, Giavalisco, 
et al. 2004a, Ouchi et al. 2004a) and Lyman a emitters 
(e.g., Hu, Cowie, & McMahon 1998, Ajiki et al. 2002, 
Kodairaetal. 2003, Rhoads et al. 2003, Hu et al. 2004). 
These surveys are sometimes not wide enough to overcome 
the effect of cosmic variance in deriving universal nature 
of these high-z populations or sometimes not deep enough 
to catch a faint signal from these distant populations. 

We carried out a large systematic survey searching for 
high-z populations with Subaru Telescope, which we call 
the Subaru Deep Field (SDF) project. Our major scien- 
tific goals of this project are to construct the largest sam- 
ples of Lyman-break galaxies at z ~ 4 — 5 and of Lyman 
a emitting galaxies at z ~ 5.7 and 6.6, and make detailed 
studies of these very high-redshift galaxy populations. We 
are also carrying out a subsequent spectroscopy and near- 
infrared imaging for these targets. In the deep optical 
imaging, we are going to take a couple of well-known 
approaches to search for these distant populations, one 
is to detect continuum break galaxies using deep multi 
broad-band imagings and the other is to isolate high-z 
populations having only the Lya emission feature by deep 
narrow-band imagings. These are the conventional meth- 
ods that have been taken in 4m class telescopes and other 
8-lOm class telescopes. However, they become very pow- 
erful with Subaru and its unique wide-field camera. 

The SDF project was approved to be one of three 
Subaru Observatory Projects. About 30 nights were al- 
located for this SDF project during 2002-2004 semesters. 
The most unique point of the SDF is that the field is 
much wider than those of other similar deep surveys. The 
Subaru prime-focus camera (Suprime-Cam; Miyazaki et 
al. 2002) has ten 2kx4k CCDs to cover a FOV as wide 
as 34' X 27'. The surveyed area is five times larger than a 
single GOODS field (Giavalisco, et al. 2004b) and even 
250 times larger than the HDF. With this very wide FOV, 
we can derive the statistical properties of distant galaxies 
such as luminosity function and angular correlation func- 
tion from large samples which are robust against cosmic 
variance. Samples of a large number of high-z galaxies en- 
able us to study different properties depending on various 
parameters, such as redshift, luminosity, and local density. 
Furthermore, the wide FOV gives us a high probability to 



discover rare objects such as the most distant galaxies. 

In the SDF project, we attempted to take the deepest 
imaging data among those having a single Suprime-Cam's 
FOV. We obtained almost 10 hours long exposure time in 
each of B, V, R, i' and z' bands. This extremely long in- 
tegration yielded large catalogues that contain more than 
100,000 objects in each band. In addition, we have also 
carried out ^ 10 hours long integrated narrow-band imag- 
ing at 8I6OA and 9210A, that correspond to dark night- 
sky windows free from OH emission lines. These are simi- 
lar wavelengths to those of other recent narrow-band deep 
surveys such as CADIS (e.g., Maier et al. 2003), LALA 
(e.g., Rhoads et al. 2003), and UH (e.g., Hu et al. 2004) 
surveys. These wavelengths corresponds to the redshifted 
Lya emission at z = 5.7 and z = 6.6, respectively. The 
initial discovery of two LAEs at z = 6.6 on the SDF was 
reported in Kodaira et al. (2003), and confirmation of an 
LAE at z — 6.33 which was detected as an i'-dropout was 
presented in Nagao et al (2004 in preparation). 

Though our survey is dedicated for high-z populations, 
we can also sample very faint galaxies in the local universe. 
For example, the SDF enables studies of galaxies down to 
M* -\- 3 around z = 1. Also our deep narrow-band images 
contain many faint emission line galaxies in foreground. 

In this paper, we present the procedures to construct 
the optical imaging data set as well as catalogues that are 
released for public. Observations are described in Section 
3, and data reduction processes are given in Section 4. We 
also show characteristics of the data such as S/N distribu- 
tion over the field, astrometric accuracies and complete- 
ness estimation (Section 6). And finally in Section 7, our 
completeness-corrected galaxy number counts are shown 
to be consistent with previous results. Throughout this 
paper, magnitudes are in the AB system. 

2. The Field 

The Subaru Deep Field (SDF) is located near 
the North Galactic Pole, being centered on 
(13''24™38.^9, -t-27°29'25."9) (J2000). We selected 
this field by imposing the following three conditions: (1) 
absence of bright stars, bright galaxies, or nearby known 
cluster of galaxies within ~ 30' from the field center, (2) 
an extreme low Galactic extinction {Ay = 0.052 at the 
center of the field; Schlegel, Flinkbeiner & Davis 1998), 
and (3) a high altitude at meridional transit from 
Mauna-Kea and thus ensuring high observing efficiency 
(the maximum altitude of the SDF is 82°). 

We conducted preliminary imaging and spectroscopic 
observations on the SDF during the commissioning runs of 
three instruments on the Subaru Telescope, Suprime-Cam 
(optical camera for the prime focus; Miyazaki et al. 2002), 
CISCO (NIR camera; Motohara et al. 2002), and FOG AS 
(optical camera and spectrograph for Cassegrain focus; 
Kashikawa et al. 2002) from 1999 to 2001. We obtained 
very deep J and K' images [Jab = 26.0, K'j^^ = 25.3; 
Maihara et al. 2001) for a 2' x 2' area near the field center. 
Some studies based on the deep NIR data were reported 
in Maihara et al. (2001); Totani et al. (2001a); Totani et 
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al. (2001b); Totani et al. (2001c); Nagashimaet al. (2002); 
Kashikawa et al. (2003). As for the Suprime-Cam im- 
ages in this prchminary survey, we obtained = 27.8 
(3hrs integration time), Vum — 27.3 (2 hrs), Rum = 27.1 
(1.5 hrs), = 26.9 (2 hrs), Z;',„ = 26.1 (1 hr), which 
are merged into the new data to construct the final im- 
ages. Some results on this wide field optical data are pre- 
sented in Ouchi et al. (2003a); Ouchi et al. (2004a); Ouchi 
et al. (2004b); Shimasaku et al. (2003); Shimasaku et 
al. (2004). 

3. Observations 

Optical imaging was made with the prime-focus cam- 
era on the Subaru Telescope, Suprime-Cam. Having 
ten 2k x4k MIT/LL CCDs, Suprime-Cam covers a con- 
tiguous area of 34' x 27' with a pixel scale of 0."202 
pixel" ^ (Miyazaki et al. 2002). The average gain of CCDs 
is 2.6e~/ADU. The widths of the gaps between adja- 
cent CCDs range over 79 — 84 pixels (or 16" — 17") in 
the north-south direction and 15 — 20 pixels (or 3" — 4") 
in the east-west direction. We imaged the SDF in five 
standard broad-band filters, B, V, R, i' , and z', and two 
narrow-band filters, NB816 (A^ = 8150A, FWHM= 120A) 
and NB921 (A^ = 9196A, FWHM= 132A) (see Figure 1). 
The two narrow-band filters, which have been carefully de- 
signed to avoid strong OH night sky emission lines, were 
used to search for Lya emitting objects at z ~ 5.7 and 
z ~ 6.6, respectively. 

The unit exposure time for each filter was determined 
so that the photon counts per pixel of the sky background 
reach an appropriate value. For the broadband filters, 
the typical unit exposure time is 900 sec for B and 240 
sec for z', being shorter for redder filters because of the 
increase in the sky brightness with wavelength. For the 
narrowband filters, the typical exposure time is 1800 sec 
for both NB816 and NB921. 

We adopted a common dithering pattern of pointings 
for all the bands. As shown in Figure 2, a full cycle of 
dithering consists of 169 pointings, which is divided into 
13 small 'circles' of 13 pointings. In the actual imaging, 
the first 13 pointings were made on a common position in 
the 13 circles (the central point in each circle). Similarly, 
the second 13 pointings were made on another common 
position in the 13 circles. The total number of pointings 
was less than 169 (a full cycle of dithering) for any filters. 
We devised the pattern of dithering, i.e., the distribution 
of 169 pointings, taking account of the following condi- 
tions (some of which are in conflict with each other). (1) 
The spread in pointings on the sky should be large enough 
to give sufficient overlaps between frames taken with ad- 
jacent CCDs. This condition is necessary, since the mo- 
saicing of frames taken with adjacent CCDs is made using 
stars common to the frames. (2) The spread in pointings 
on the sky should be small, in order to reduce areas of low 
S/N ratios near the edges of the mosaiccd image. (3) The 
mosaiced image should have sufficiently uniform S/N ra- 
tios. (4) In the course of dithering, no sky position should 
fall onto the same dead-pixel area more than once. Each 



CCD has a number of long (but narrow) lines of dead 
pixels. To meet condition (4), we adjusted the dithering 
positions in the horizontal {X) and vertical {Y) axes so 
that there are no duplicated X values or Y values in the 
169 positions. 

As described below, the final images and catalogs are 
based only on the data with a seeing size of < 0."98. The 
fraction of the frames satisfying this seeing condition at- 
tains at least 70% in exposure time for any of the filters. 
The total exposure times, limiting magnitudes {3a on a 2" 
aperture), and observing dates are summarized in Table 1. 
Table 1 presents only the data good enough to be used 
to make final, mosaiced images. Low-quality data taken 
either in poor seeing (worse than 0."98) or in very low 
transparency arc not listed in the table. All the raw data, 
including those with low quality, are available at SMOKA 
(Subaru Mitaka Okayama Kiso Archive) Science Archive . 

In constructing the final mosaiced images in the B, V, 
i' , and z' bandpasses, we have co-added the data with 
1 — 3 hrs integration times taken during the commission- 
ing observations of Suprime-Cam in 2001. The i'-band 
data were taken separately over several observing runs 
because a high-z supernova search project (Yasuda et al. 
in preparation) in this field requires intermittent observa- 
tions in this band. This repeated observation also enables 
the search for time- variable objects. 

In the B, V, and R observations, one of the three 
standard photometric standard fields, SA104, SA107, and 
SAllO (Landolt 1992), was observed for flux calibra- 
tion, when the nights were photometric. Similarly, in 
the i', z', NB816, and NB921 imaging, at least one of 
the five spectrophotometric standard stars, HZ-21, HZ- 
44, GD153, P177D, and P330E, was observed. 

4. Data Reduction 

We reduce all the data with good quality using a 
pipeHne software package (Ouchi et al. 2003b) whose 
core programs were taken from IRAF and from soft- 
ware for mosaic-CCD data reduction developed by Yagi 
et al. (2002). The package includes bias subtraction, flat- 
fielding, correction for the image distortion due to the 
prime-focus optics, PSF matching (by Gaussian smooth- 
ing), sky subtraction, and mosaicing. Bias subtraction 
and flat fielding were processed in the same manner as 
for the conventional single chip CCD. The flat frame was 
made from the median frame of a number of normalized 
object frames. The distortion correction was performed 
with 4-th order polynomial transformation (Miyazaki et 
al. 2002) with ~ 0.5 pix rms residuals. The PSF matching 
was made for all frames so that the final mosaiced im- 
ages in the seven bandpasses have a common PSF FWHM 
value of 0."98 through smoothing the images with ade- 
quate Gaussian kernels. The global sky background was 
determined from bilinear interpolation of the typical sky 
counts of 64 x 64 pix (12". 8 x 12". 8) meshes. In mosaic- 
ing, first, the relative positions and relative throughputs 

* http://smoka.nao.ac.jp/ 
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between frames taken with different CCDs and at dif- 
ferent exposures are calculated using stars common to 
adjacent frames. Then, all the frames arc merged into 
a single, large image by adjusting relative positions and 
stacking with weights according to the relative through- 
put of frames. The geometry differences among seven 
band images were corrected with IRAF's geomap and 
geotrcin tasks. About 1000 objects were selected evenly 
over the image using SExtractor, and then spatial trans- 
formation law was determined from the iteration of 3rd 
or higher order polynomial fitting with 3a clipping un- 
less the rms residuals were smaller than O.lpix. The final 
mosaiced images for the seven bandpass are available at: 
http://soaps.naoj.org/sdf/. The size of the final images 
is 8900 X 11000 pixels (corresponding to 30' x 37') after 
removal of very low-S/N regions near the edges of the im- 
ages. In the future, we plan to make another set of SDF 
image data which consist of the images with good seeing 
sizes only in spite of shorter net integration times. 

Object Detection and Photometry 

We perform object detection and photometry by run- 
ning SExtractor version 2.1.6 (Berlin & Arnouts 1996) on 
the images. The object detection is made for all seven 
images (five broadbands and two narrowbands) indepen- 
dently. For all objects detected in a given bandpass, their 
magnitudes and several parameters in the other six band- 
passes are also measured at exactly the same positions 
as in the detection-band image, using the 'double image 
mode' of SExtractor. We detected objects that have 5 con- 
nected pixels above 2a of the sky background rms noise 
and made photometric measurements at 2a level in all the 
catalogs. Overlapping objects are deblended into distinct 
objects if the contrast between integrated intensity asso- 
ciated with each object is > 0.005. Aperture photometry 
was done with 2"4) and 3"0 apertures. Our input param- 
eters on running SExtracter and the full list of the output 
parameters contained in the catalogs can also be browsed 
at: http://soaps.naoj.org/sdf/ 

4-. 2. Photometric Calibration 

We calculate the photometric zero points of the images 
from the photometry of standard stars with a lO"^ aper- 
ture, taking account of the differences in airmass between 
the standard-star frames and the fiducial SDF frames. We 
check the photometric zero points using the stellar col- 
ors calculated from the spcctrophotomctric atlas of 175 
Galactic stars given in Gunn & Strykcr (1983). Since the 
FOV of Suprime-Cam is large, more than 100 bright (but 
not saturated) stars with i' ^ 23 are commonly detected 
in mosaiced images of the seven bands. We compare the 
location of the stellar sequence in various two-color planes 
between our data and Gunn & Striker's atlas, finding that 
the colors of stars in our data are in agreement with those 
of Gunn & Striker's stars within ~ 0.05 mag. Small offsets 
up to ~ 0.05 mag seen in some combinations of bandpasses 
may indicate that the nights were not strictly photometric 
when the standard stars were taken for those bandpasses. 
We correct the zero points of all the bandpasses by adding 



a small constant (at most ~ ±0.05 mag) to the original 
value, so that the observed colors of stellar objects match 
completely the synthetic colors of Gunn & Stryker's stars. 
In any case, the errors in the photometric-zero points of 
our final images would be less than 0.05 mag. 

4-3. Limiting Magnitudes 

We measure the 3a limiting magnitudes of the images, 
which are defined as the 3a of the sky noise on a 2"(j> diam- 
eter. For each image, we measure sky counts in a number 
of 2"(j) apertures which are located at randomly selected 
positions over the image. Then, we draw a histogram of 
the sky counts, and fit a Gaussian function to the his- 
togram to obtain the la noise. When we fit a Gaussian 
function, we do not use the positive tail of the histogram 
which is affected by objects. The limiting magnitudes thus 
estimated are presented in Table 1. Note that we did not 
simply calculate the 3a limiting magnitudes by scaling the 
la noise in one pixel to a 2"-diameter area, since such a 
simple scaling gives a very optimistic limiting magnitude 
as discussed in Furusawa (2002) and Labbe et al. (2003). 

4.4- Astrometry 

The coordinates of the final, mosaiced images agree very 
well with the tangential coordinates, i.e., coordinates ob- 
tained by projecting the sky onto a plane. Our mosaicing 
procedure using a large number of stars common to adja- 
cent unit frames ensures very small random errors in the 
relative position (and relative throughput) over the entire 
mosaiced images; the relative errors are typically less than 
1 pixel (or 0."2). 

We then make absolute astrometric calibration as fol- 
lows. We calibrate the coordinates (World Coordinate 
System) of the images using the objects common to cat- 
alogs based on our data and the 2MASS All-Sky Catalog 
of Point Sources (2MASS-PSC) (Cutri et al. 2003). 
Calibrations are done for each band separately. The po- 
sitional accuracies of 2MASS-PSC sources are as good as 
70-80mas over most of the magnitude range and 200mas 
even at the faintest magnitudes of K = 16. We use 385 
2MASS-PSC sources which are located in the SDF area. 
They arc uniformly distributed over the image. On the 
other hand, we have selected un-saturated SDF objects 
which have over 50 connected pixels above lOcr of the sky 
background. Then we took cross-correlation between the 
positions of SDF objects and 2MASS-PSC sources. We 
find 277-290 objects are matched between two catalogs 
after rejection of some obvious disagreements due to ob- 
ject blending. The faintest magnitudes of matched ob- 
jects are 23.03, 21.85, 20.88, 19.89, 18.75, 18.89 and 18.66 
for B, V, R, i', z', NB816 and NB921 band, respec- 
tively. Calculations of WCS have been done using IRAF's 
images . imcoords . ccmap task. 

The fitting procedure made a small residuals as shown 
in Figure 3. The absolute position accuracy was found 
to be 0.224, 0.258, 0.266, 0.259, 0.214, 0.222 and 0.209 
arcsec, and rms errors arc 0.154, 0.198, 0.205, 0.196, 0.152, 
0.160 and 0.144 arcsec in right ascension, 0.163, 0.167, 
0.171, 0.171, 0.150, 0.154 and 0.152 arcsec in declination 
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for B, V, R, i', z', NB816 and NB921 band, respectively. 
Figure 3 also shows that there is no systematic dependence 
on the magnitude and no significant offsets either in right 
ascension and declination. 

5. Catalogs 

We make seven catalog set for each detection bandpass. 
For instance, the B-band detected catalog set contains 
seven measured catalogs which list the position, magni- 
tudes, and several parameters concerning size and shape, 
of objects measured in the i?-band image, as well as the 
magnitudes in the other six bandpasses measured for the 
same position and aperture as in the S-band image. 

The efficiency and reliability of object detection and 
photometry are significantly lower than the average 
around very bright stars (due to bright haloes and satura- 
tion trails) and near the edges of the images (due to low 
S/N ratios). We carefully defined these low-quality re- 
gions ( "masked regions" ) which correspond to low-quality 
regions such as frame edges and the area around very 
bright stars with saturation trails. The final catalogs are 
not masked, though it is better for most of cases to remove 
objects falling onto these regions. Figure 4 plots the S/N 
distribution map on the final image in the case of the i'- 
band. The sky noise values are evaluated in the same man- 
ner as was described in Section 4.3 for each of 355 x 355 
pixel meshes with steps overlapping a half size of the mesh. 
Then the S/N values are calculated in the case of the ob- 
jects with the limiting magnitude of i' = 27.43. Darker 
regions imply highcr-5/A'^ regions. Figure 4 demonstrates 
that our dithering pattern gives a fairly uniform distri- 
bution of S/N ratios over the region. The almost ax- 
isymmetric distribution rcfiects the vignetting pattern of 
prime-focus optics and the slightly lower S/N region near 
the edge in the left-hand side would be due to the auto- 
guider probe which sometimes shadowed the FOV. The 
lower left-hand corner with considerably lower S/N ratios 
refiects the low quantum efficiency of the CCD which cov- 
ered this region; the quantum efficiency of this CCD is 
about 60% those of the others due to an incomplete AR 
coating (Miyazaki et al. 2002). Similar S/N maps have 
been obtained for the other filters. 

The number of objects contained in the final catalogs 
are 214,733 {B), 198,466 (V), 209,477 (R), 198,581 (i'), 
187,647 (z'), 195,369 (NB816), and 197,363 (NB921). 
If we remove objects that are positiond in the masked 
regions, the numbers of objects arc reduced to 162,611 
{B), 150,991 {V), 160,129 (E), 153,030 (i'), 126,332 (z'), 
133,242 (NB816), and 119,493 (NB921). The catalogs 
as well as the file which defines the masked regions are 
available at: http://soaps.naoj.org/sdf/. 

6. Completeness 

In this section, we describe the simulation carried out 
to estimate the completeness of our catalogs. For that 
purpose, we distribute artificial objects on the original im- 
ages. To characterize an ensemble of the artificial galaxies, 



we have to assume the distributions of magnitude, surface 
brightness, luminosity profile. These parameters strongly 
influence on the estimate of sample completeness. Though 
the completeness of galaxy sample should be determined 
as a function of both magnitude and surface brightness, 
our sample is so deep that it is not sure that the assump- 
tions for these faint galaxies are valid. Therefore, we avoid 
an inaccurate estimate for extended objects whose com- 
pleteness significantly depends on surface brightness, and 
alternatively, we present more secure evaluation of the 
sample completeness for point sources. The extended ob- 
jects would have brighter complete limits than those we 
show here. 

We used IRAF tasks of star list and mkobjects to 
make artificial stars on the original images. Artificial 
stars have the same FWHM (0.98arcsec) as the real im- 
ages, and are randomly distributed on the original images. 
This random layout sometimes places artificial objects on 
real objects, which causes a heavy blending and incorrect 
photometry. Thus, we neglect an artificial object if the 
distance from it to the nearest real object is less than 
three times the FWHM of the real object. This rejection 
procedure means that we evaluate the object detection ef- 
ficiency on pure sky background regions, and thus may 
give slightly optimistic estimation for completeness be- 
cause object blending actually reduces the efficiency. We 
extract artificial stellar objects with SExtractor using the 
same parameter set as adopted in the catalog construction 
(Section 5). Then we measure the fraction in number of 
detected objects to the input sources. We generate 50,000 
artificial objects for a realization and repeat it four times 
to obtain a statistically robust number (~ 10,000) of ob- 
jects for each 0.5 magnitude bin. Figure 5 presents the 
resulting completeness estimates for each broadband as a 
function of total magnitude (MAGJVUTO magnitudes from 
SExtractor). 

7. Galaxy Counts 

In this section, we compare the galaxy number counts 
in the SDF with those given in the literature for each 
broad band. When deriving the galaxy counts from our 
catalogs, we remove bright (< 24 mag) Galactic stars from 
the catalogs, by imposing the following criteria: 

/star > 0.99 (mag < 20), (1) 
log(FWHM) < (24 - mag)/ 150(20 < mag < 24). (2) 

where /star is the stellarity index, an output from 
SExtractor, FWHM is in arcscc, and mag is the total 
magnitude. We do not make star/galaxy separation for 
objects fainter than 24 mag, since in those magnitudes 
a significant fraction of galaxies are too small to securely 
discriminate from stars under the seeing size of our images 
(0."98). However, the surface density of stars is expected 
to be much lower than that of galaxies in those faint mag- 
nitudes. We then correct the raw galaxy counts for detec- 
tion completeness. On the other hand, our galaxy catalogs 
can include false spurious objects. We estimate the con- 
tamination by spurious objects by running SExtractor on 
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negative images in the same manner as for the original 
images. We find that the contamination is as low as 2% 
even in the faintest magnitude bins of our number counts. 
Thus, we do not correct our counts for this small contam- 
ination by spurious objects. 

Figure 6(a-e) plots the galaxy counts derived from our 
data together with those taken from the literature. The 
filled circles and open circles correspond, respectively, to 
the raw counts and the corrected counts. Since the raw 
counts include Galactic stars, they are higher than the cor- 
rected counts at the bright end of the plots. These counts 
are also presented in Table 2. The solid lines shown at the 
bright end for the B, V, and R plots are predicted star 
counts toward the SDF based on the stellar population 
synthesis model for the Galaxy by Robin et al. (2003). 
We have confirmed that the star counts calculated from 
our data (< 24 mag) are consistent with the predicted 
counts in all three bandpasses. Note that i'- and z'-band 
predictions are not available. 

The galaxy number counts in the SDF data are reliable 
up to ~ 28 mag for B, ~ 27.5 mag for V, R, and i' , and 
~ 26.5 mag for z' , beyond which the completeness drops 
below 30%. At these faint-end magnitudes, the counts 
reach ~ 3 x 10^ galaxies per 0.5 mag per square degree for 
all five bandpasses. Comparison between our counts and 
those given in the literature shows broad agreement over 
the entire overlapping magnitude range for all bandpasses. 
The SDF galaxy counts cover as wide as 8 magnitudes 
from 20 mag up to ~ 28 mag. 

8. Summary 

The paper describes an overview of the SDF project and 
the procedures to complete the optical imaging data and 
catalogues. We have carried out a deep optical imaging in 
seven filters, B, V, R, i', z' , NB816, and NB921 over the 
34' X 27' in the SDF. The limiting magnitudes measured 
at 3cr on a 2" aperture are B = 28.45, V = 27.74, R = 
27.80, i' = 27 A3, z' = 26.62, NB816 = 26.63, and NB921 = 
26.54 in the AB system. Object catalogs constructed from 
these extremely deep data contain more than 10^ objects 
in each bandpass. The galaxy number counts corrected 
for detection completeness are found to be consistent with 
previous results in the literature. 

Scientific results based on the data will be shown in 
the forthcoming papers. We hope that the public released 
SDF data will help in advances of our understanding for 
the high-z universe as well as for a variety of astronomical 
topics. 

We deeply appreciate the devoted technical and man- 
agemental support of the Subaru Telescope staff to this 
long-term project. The observing time for this project 
was committed to all the Subaru Telescope builders. We 
thank the referee for helpful comments that improved 
the manuscript. The research is supported by the Japan 
Society for the Promotion of Science through Grant-in- 
Aid for Scientific Research 16740118. 



No. ] SDF Imaging Data 

Table 1. Summary of Optical Imaging Data for SDF. 
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band 


cxp.time 


seeing size 


miim* 






Date of Observations 




(min.) 


(arcsec) 


(mag) 


(mag/count) 






B 


595 


0.98 


28.45 


34.780 


162,611 


2003, Mar. 31, Apr. 2 


V 


340 


0.98 


27.74 


34.584 


150,991 


2003, Apr. 30 


R 


600 


0.98 


27.80 


34.315 


160,129 


2003, Mar. 30, Apr. 1/2/25/29/30 


i' 


801 


0.98 


27.43 


34.253 


153,030 


2002, Apr. 11/14, May 6, 


z' 


504 


0.98 


26.62 


33.001 


126,332 


2003, Mar. 31, Apr. 2/24/25/29/30 
2002, Apr. 9/14, 2003, Mar. 7, Apr. 1/28 


NB816 


600 


0.98 


26.63 


32.880 


133,242 


2003, Apr. 28/29/30 


NB921 


899 


0.98 


26.54 


32.520 


119,493 


2002, Apr. 9/11/14, May 6, 

2003, Mar. 7/8, Apr. 24 



* The 3ct limiting magnitudes in the AB system within 2" diameter aperture. 
^ The magnitude zero point. 

t The number of detected objects out of the masked regions. 



Table 2. Number counts on SDF.^* 



Magnitude 




B 






V 






R 






i 








2 


./ 




18, 


,250 


2, 


.076 


1, 


.882 


2, 


,246 


2, 


,169 


2, 


.350 


2, 


.169 


2, 


,607 


2, 


.280 


2, 


,642 


2. 


,234 


18, 


,750 


2, 


.155 


1, 


.155 


2, 


,386 


2, 


,076 


2, 


.597 


2, 


,269 


2, 


,739 


2 


.394 


2, 


,802 


2. 


,523 


19, 


,250 


2, 


,394 


1, 


.792 


2, 


,597 


2, 


,222 


2, 


.695 


2, 


,291 


2, 


,848 


2 


.504 


2, 


,885 


2. 


,592 


19, 


,750 


2, 


,368 


1, 


.909 


2, 


,674 


2, 


,301 


2, 


.868 


2, 


,597 


2, 


,955 


2, 


.678 


3, 


,100 


2. 


,919 


20, 


,250 


2, 


.477 


2, 


.155 


2, 


,805 


2, 


,570 


2, 


.990 


2, 


.812 


3, 


,149 


3 


.000 


3, 


,282 


3. 


,185 


20, 


,750 


2, 


,682 


2, 


.394 


3, 


,015 


2, 


,848 


3 


.175 


3, 


,055 


3, 


,360 


3 


.258 


3, 


,486 


3. 


,399 


21, 


,250 


2, 


,848 


2, 


.651 


3, 


,124 


3, 


,019 


3 


.361 


3, 


,279 


3, 


,529 


3 


.472 


3, 


,644 


3. 


,595 


21, 


,750 


3, 


,056 


2 


.941 


3, 


,322 


3, 


,233 


3 


.544 


3, 


,477 


3, 


,697 


3 


.651 


3, 


,847 


3. 


,808 


22, 


,250 


3, 


,265 


3, 


.195 


3, 


,537 


3, 


,479 


3 


.743 


3, 


,698 


3, 


,882 


3 


.848 


3, 


,959 


3. 


,927 


22, 


,750 


3, 


,535 


3, 


.480 


3, 


.747 


3, 


,711 


3 


.918 


3, 


,891 


4, 


,011 


3 


.986 


4, 


,113 


4. 


,101 


23, 


,250 


3, 


,774 


3 


.747 


3, 


,976 


3, 


,955 


4, 


.087 


4, 


,068 


4, 


,177 


4, 


.181 


4, 


,261 


4. 


,261 


23, 


,750 


4, 


.048 


4, 


.034 


4, 


,194 


4, 


,192 


4, 


.290 


4, 


,293 


4, 


,354 


4, 


.367 


4, 


,417 


4. 


,430 


24, 


.250 


4, 


.285 


4, 


.299 


4, 


,410 


4, 


,429 


4, 


.480 


4, 


.498 


4, 


,516 


4, 


.549 


4, 


,559 


4. 


,599 


24, 


,750 


4, 


,483 


4, 


.501 


4, 


,590 


4, 


,617 


4, 


.649 


4, 


,674 


4, 


,661 


4, 


.698 


4, 


,694 


4. 


,785 


25, 


,250 


4, 


.652 


4, 


.677 


4, 


,739 


4, 


,780 


4, 


.770 


4, 


,807 


4, 


,781 


4, 


.846 


4, 


,790 


4. 


,972 


25, 


,750 


4, 


.779 


4, 


.819 


4, 


,846 


4, 


,942 


4, 


.871 


4, 


,942 


4, 


,873 


5, 


.007 


4, 


,863 


5. 


,171 


26, 


,250 


4, 


,882 


4, 


.961 


4, 


.932 


5, 


,133 


4, 


.940 


5, 


,089 


4, 


,932 


5, 


.179 


4, 


,885 


5. 


,332 


26, 


,750 


4, 


,952 


5, 


.119 


4, 


.973 


5, 


,300 


4, 


.981 


5, 


,246 


4, 


,942 


5, 


.329 










27, 


,250 


4, 


,989 


5, 


.287 


4, 


,953 


5, 


,418 


4, 


.966 


5, 


,357 


















27, 


,750 


4, 


,973 


5, 


.404 



































* The first columns show the raw counts in log(/deg'^/0.5mag), and the second column 
show the counts corrected for completeness and star counts for each band. 
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SDF Imaging Data 
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Fig. 1. The band response curves of SDF optical imaging. 
These are convolved with the CCD sensitivities, instrument 
and atmospheric transmissions. Thin black line shows the OH 
night sky line spectrum. 
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Fig. 2. The dithering pattern of observing pointings of SDF. 
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Fig. 3. Astrometry accuracy of SDF in the case of i'-band. 
Aa and AS are calibration residuals between SDF objects and 
2MASS-PSC sources in the directions of right ascension and 
declination, respectively. 
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Fig. 4. The S/N map of SDF in the case of j'-band. Thick 
(red) line shows the trimmed area defined in the "masked 
region" file. 



SDF Imaging Data 




Fig. 5. The complctess for our final SDF sample as a fune- 
tion of magnitudes. 
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Fig. 6. The galaxy number counts in B-band. Open (red) 
circles show the raw number counts from the SDF catalog, 
filled (red) circles show the galaxy number counts corrected 
with completeness and subtracted star counts, and other sym- 
bols from the literature. Solid line shows the predicted star 
counts of Robin et al. (2003). The error bars shown in the 
SDF counts are based on the simple Poisson errors. 
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Fig. 7. Same as in Figure 6 but in V-band. 
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Fig. 8. Same as in Figure 6 but in R-band. 
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Fig. 9. Same as in Figure 6 but in i'-band. 
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Fig. 10. Same as in Figure 6 but in z'-band. 



